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Abstract 

The search for quartic anomalous gauge couplings at LEP requires appropriate 
predictions for the radiative processes e + e~ — > 1/P77, e + e~ — > qqj'y and e + e~ — > 
4 fermions+7. Matrix elements are exactly computed at the tree level, and the 
effects of anomalous couplings and initial-state radiation are included. Comparisons 
with results and approximations existing in the literature are shown and commented. 
Improved versions of the event generators NUNUGPV and WRAP are made available for 
experimental analysis. 
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1 Introduction 



Despite of the striking success of the Standard Model (SM) in accommodating the pre- 
cision data collected at high-energy colliders, important tests of the theory, such as the 
non-abelian nature of the gauge symmetry and the mechanism of electroweak symme- 
try breaking, are still at a beginning stage. To this end, gauge-boson self interactions 
play a key role. Presently, triple gauge couplings are being probed at LEP [1] and Teva- 
tron [2], while only very recently direct measurements of quartic couplings became avail- 
able through the study of radiative events at LEP [3]- [6]. Actually, events with one or 
two isolated, hard photons are analyzed at LEP to search for anomalies in the sector of 
quartic gauge-boson couplings. Only vertices involving at least one photon can be con- 
strained since quadrilinear interactions containing four massive gauge bosons give rise to 
a three massive gauge boson final state and are therefore beyond the potentials of LEP 
due to lack of phase space. The processes considered in the experimental analyses are 
e + e~ — > W + W~y ) e + e~ — > Z77 and e + e _ — > z/z/77 [3]- [6]. The W + W~ / y signature, which 
yields a four fermion plus gamma final state, is interesting in order to test WW'-ff and 
WW Zy vertices. For the Z77 events, the final states due to the hadronic decays of the Z 
boson with two jets and two visible photons are selected to probe the purely anomalous 
vertex ZZ77, which is of particular interest being absent in the SM at tree level. The final 
state with two neutrinos and two acoplanar photons allows to study the quartic WW^yy 
and ZZ'-f-f interactions. 

For tests of the SM and searches for new physics beyond it, anomalous quartic gauge- 
boson couplings are important because, as widely discussed in the literature [7], they 
offer a window on the mechanism of symmetry breaking and it is possible to imagine 
extensions of the SM that alter the quartic vertices without modifying the trilinear in- 
teractions. For this reason, phenomenological analyses of anomalous quartic interactions 
have already appeared in the literature, by considering the potentials of 77 [8], ey [9] 
and e + e" [10]- [13] colliders and very recently also of hadron machines [14]. The tight- 
est constraints on quartic anomalous couplings come from electroweak precision data, as 
discussed in ref. [14]. However, present measurements at LEP2 can investigate photonic 
quartic couplings for the first time in a direct way through the study of the final states 
above discussed. In the light of these experimental analyses, the aim of this letter is to 
present exact SM calculations of the processes e + e~ — > 4 fermions (4/) + 7, e + e _ — > qqyy 
and e + e" — > z/z/77, including the effects of quartic anomalous gauge couplings (QAGC) 
and of the most important radiative corrections due to initial-state radiation (ISR). The 
sensitivity of LEP searches is studied in comparison with approximate results existing in 
the literature. It turns out that the Z77 approximation works well for the (7(777 and z/z/77 
channels when appropriate cuts around the Z mass are imposed, while the WWy approx- 
imation can significantly differ from the exact 4/ + 7 calculation, even in the presence 
of invariant mass cuts around the W mass. It is also shown that the z/z/77 fi na l state 
with appropriate cuts on the recoil mass can be successfully exploited to extract limits 
on neutral QAGC, as a complementary channel to the (7(777 process. 
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2 Theoretical approach 



The theoretical framework of interest is the formalism of electroweak chiral lagrangians. 
In such a scenario, QAGC involving four massive gauge bosons emerge as operators of 
dimension four at next-to-leading order, while QAGC with at least one photon originate 
from six (or higher) dimensional operators at next-to-next-to-leading order. They are said 
genuinely anomalous if they do not induce new trilinear gauge interactions. 

Anomalous WW'yy and vertices were originally introduced in ref [8]. In this 

paper, the authors show that, by assuming C and P conservation and further imposing 
U(l)em g au g e invariance and SU(2) C custodial symmetry, two independent Lorentz struc- 
tures contribute to WW^yy and interactions according to the following lagrangians 

Co = -^F, U F^W«-W a 

£c = — ^V^w^-Wp, (1) 

where F^ u = d {l A v — d^A^ is the electromagnetic field tensor, and W is a SU (2) triplet 
describing the W and Z physical fields, i.e. 



W 



\ ZJ cos 9 W J 



cos9 w being the cosine of the weak mixing angle. In eq. (1) a and a c are (dimensionless) 
anomalous couplings, divided by an energy scale A, which has the meaning of scale of new 
physics. Generally speaking, A is in principle unknown and model-dependent. However, 
the ratios a«/A 2 entering the phenomenological lagrangians can be meaningfully extracted 
from the data in a model-independent way. 

An anomalous WW Z^ interaction has also been advocated in refs. [9]-[ll] . The cor- 
responding phenomenological lagrangian reads as follows 



e 2 a r 



£n = -^ 2 e ljk W; a WlW ka F^ , (2) 

where W^ v is the SU(2) field strength tensor, and a n /A 2 the anomalous coupling. Al- 
though U(l) em and SU{2) C conserving, the lagrangian of eq. (2) violates C and CP, as 
noticed in the literature [13]. The anomalous couplings a ,a c ,a n entering eqs. (l)-(2) 
are presently constrained by LEP collaborations [3]-[6]. Under the assumptions of lo- 
cal U(l) em invariance and global custodial SU(2) C symmetry, two additional operator 
structures which violate P have been very recently proposed in ref. [15]. 

A complete and general analysis of photonic quartic couplings has been performed in 
ref. [13]. In this paper it is shown that, by imposing C, P and U(l) em invariance, nine 
independent Lorentz structures do contribute to WW'yy, ZZ77 and WWZ^y vertices. It 
is further demonstrated by the authors that, by embedding these structures in SU(2) x 
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U(l) gauge invariant and SU(2) C symmetric combinations, fourteen C and P conserving 
operators are allowed, with kj parameters, which parameterize the strength of anomalous 
couplings. 

Following ref. [13], the nine independent operator structures contributing to the ver- 
tices analyzed at LEP have been implemented directly at the lagrangian level in the ALPHA 
algorithm [16], according to the formula 

J^qagc = 

W 1 + W 2 + Z 1 + z 2 

W z + W z + W z + W z + W z , (3) 
where the Lorentz structure of the operators is given by: 



W 1 = 


a^F^W+W-o 




w 2 = 


a w2 F^F^W + »W- 


+ h.c. 


Zi = 


a zl F^F^Z p Z<> 




z 2 = 


a z2 F pu F^Z»Z p 




W Z = 


a wz0 F, u Z^W;W-t 


> 


w z = 


a wzc F^Z^W + ^W- 


+ h.c. 


w z = 


a^F^W+^Z^W- 


" +h.c. 


w z = 


a wz2 F pv W^Z^W; 


" +h.c. 


w z = 


a wz3 F^W + ^Z p W- 


" + h.c 



(4) 



the cii being coefficients of dimension M~ 2 . 

The first two structures refer to WW'yy interactions, the third and fourth ones to 
ZZ^ff interactions, while the remaining five affect the WW Z^ vertex. It is worth noticing 
that, by virtue of eq. (4), both parameterizations available in the literature for QAGC, 
namely the parameterization in terms of ao, a c , a n couplings and the one in terms of kj 
coefficients, can be obtained by means of appropriate relations between the a, parameters. 
For example, assigned a , a c and a n , the following relations for the coefficients of eq. (4) 
hold 

e 2 

&wl 8A 2 " ^0 
e 2 

azl ~ 16cos 2 u ,A 2 a ° 
e 2 

a ^ = "lecos^^A 2 ° c ( 5 ) 
_ ' 

awzc = 1 16 cos 6»„A 2 ° n 
_ ' 

awz2 ~ % l&cose w K 2 ttn 

_ • (p 1 
&wzZ — —1 16cos 8wA 2 a„ . 

The implementation in ALPHA has been carefully cross-checked by an analytical cal- 
culation of the scattering processes WW — > 77, ZZ — > 77 and WW — > Z7 (and its 
permutations), finding perfect agreement. 

In order to provide valuable tools for experimental data analysis, the version of ALPHA 
including the lagrangian of eq. (4) has been interfaced to the Monte Carlo generators 
NUNUGPV [17] and WRAP [18], for a complete study of QAGC in radiative events at LEP. 
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By NUNUGPV and WRAP 
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Figure 1: The cross sections in the absence (solid line) and in the presence (dashed and dotted 
lines) of anomalous parameters ao and a c (ao/A 2 = a c /A 2 = 0.05 GeV -2 ) for the processes 
e + e~ — > Z77, WWj, 9577, 1^77, ud^v^, as functions of the cm. energy in the LEP2 energy 
range. 



The z/z/77 final state can be simulated by means of NUNUGPV, which includes exact SM 
matrix elements and p r dependent QED Structure Functions (SF) to account for ISR. The 
other processes of interest at LEP can be studied by using WRAP, which is based on the 
exact matrix elements for the (charged-current) reactions e + e~ — > 4/+7 and the process 
e + e~ — > qq r y r ) and the same treatment as NUNUGPV for ISR. By means of WRAP, predictions 
for the inclusive final states WW^ and Z77 can be also obtained, especially in order to 
compare with results existing in the literature treating the W and Z bosons in the on-shell 
approximation, as those of refs. [11, 13]. 

3 Numerical results and discussion 

In order to test the implementation of QAGC, comparisons with the available results have 
been performed. A first comparison has been done with the plots published in ref. [13] 
for the dependence of the WW^j and Z77 cross sections on the anomalous parameters 
k{ at y/s = 200 GeV, using the same set of input parameters and cuts. The results of 
such a comparison between the results of ref. [13] and the predictions of the Monte Carlo 
WRAP show a very satisfactory agreement, testifying a correct implementation of WW'yy, 
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Figure 2: The sensitivity of the qqj^y and 1/1/77 final states to a>o and a c anomalous couplings 
in comparison with the Z77 approximation, at \/i =200 GeV. 



ZZ77 and WW Z^i couplings in the present calculation. 

As far as the parameterization in terms of a and a c anomalous couplings is concerned, 
two further comparisons have been performed. The first one has been done with the 
results of ref. [13] for the dependence of the WW^ cross section on ao, finding perfect 
agreement. The second one has been performed with the results published in ref. [11] for 
the dependence of the WWj and Z77 cross sections on a and a c parameters, using the 
same set of cuts. The predictions obtained with WRAP differ from those of ref. [11]. More 
precisely, an opposite sign is present for the relative effects of a and a c on the WW^y 
cross section, as noticed in ref. [15] and confirmed [19] by the authors of ref. [11], whereas 
this is not the case for the Z77 final state. Work is in progress in order to clarify this 
point [19]. 

For the results shown in the following, the input parameters used are: 

G F = 1.16637- 1CT 5 GeV" 2 M z = 91.1867 GeV 

M w = 80.35 GeV sin 2 6 W = 1 - /Mf (6) 

T z = 2 GeV T w = 2.04277 GeV 

The cuts adopted are: E 1 > 5 GeV, 15° < -# 7 < 165° as detection criteria for the 
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a D /A 2 = -0.1 GeV~ 2 
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Figure 3: Distribution of the energy of the most energetic photon (-E71), the energy of the 
second most energetic photon (-E72)) the recoil mass (M reco ii) and the polar angle of the most 
forward photon (max(| cosi?i|, | cos#2|))) for the uujy final state, with ao/A 2 = —0.1 GeV -2 
and =200 GeV. Dotted histograms: SM calculation; solid histograms: full Z/P77 calculation 
for ao/A 2 = —0.1 GeV -2 ; dashed histograms: ^^77 prediction for ao/A 2 = —0.1 GeV -2 . 



observed photons, together with a photon-final state charged fermion separation cut of 
5°. Invariant mass cuts of the kind 80 GeV < M q q < 100 GeV for the e + e - — > qqyy 
process and M M j >10 GeV for e + e - — > ud/iP^ process are also imposed. 

Figure 1 shows the total cross sections, as obtained by means of the event generators 
NUNUGPV and WRAP, of all the radiative processes accessible to the LEP investigation, as 
functions of the center of mass (cm.) energy in the LEP2 energy range. The solid lines 
correspond to the pure SM predictions; the dashed and dotted lines illustrate the effects 
of anomalous couplings a /A 2 = 0.05 GeV -2 and a c /A 2 = 0.05 GeV -2 , respectively. The 
unitarity violations induced by the anomalous couplings are clearly visible, their effects 
growing as the cm. energy increases. 

In Fig. 2 the sensitivity of the final states qq'yy (dotted line) and z/z/77 (dashed line), 
as obtained by means of the full calculation, to the and a c parameters is shown at the 
cm. energy of 200 GeV, in comparison with the Z77 approximation (solid line). The 
predictions for the qq'yy channel refer to a cut on the invariant mass of the jet-jet system 
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a c /A 2 = -0.1 GeV~ 2 
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Figure 4: The same as Fig. 3, for a c /A 2 = -0.1 GeV~ 2 . 



around the Z mass (80 GeV < M q q < 100 GeV), while for the z/z/77 final state a cut 
on the recoil mass again around the Z mass (80 GeV < M reco n < 120 GeV) has been 
imposed. These cuts, also typically adopted in the experimental analysis, are required in 
order to compare with the Z77 approximation, which, according to the results available 
in the literature [11, 13], assumes the Z boson as an on-shell particle. Therefore, this 
comparison allows to quantify the effects due to the 7-Z interference in the qq'-fj channel 
and to W-Z interference in the z/z/77 one ) as wen as the effects of the off-shellness of the 
Z boson, in the extraction of limits on the QAGC. It can be seen from Fig. 2 that the 
Z77 calculation is a very good approximation of the full prediction, which includes in- 
terference contributions. A detailed numerical investigation shows an agreement between 
the integrated cross sections at the per cent level as a function of ao and Oc variations 
inside the presently allowed experimental constraints, thus illustrating the reliability of 
the Z77 approximation in view of the expected experimental precision. 

This conclusion also holds for the differential distributions mostly sensitive to QAGC 
and considered in the experimental studies, as can be noticed from Figs. 3-4, showing 
the energy of the most energetic photon (E 7 i), the energy of the second most energetic 
photon (E^) , the recoil mass (M rcco ii) and the polar angle of the most forward photon 
(max(|cos#i|,|cos# 2 |)), for a /A 2 = -0.1 GeV" 2 (Fig. 3) and a c /A 2 = -0.1 GeV" 2 
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Figure 5: The sensitivity of the ud\iv^ final state to ao, a c and a n anomalous couplings in 
comparison with the WWj approximation, at y/s =200 GeV. 



(Fig. 4). In Figs. 3-4 the predictions obtained by means of the complete calculation for 
the z/z/77 final state, including W-Z interference, (solid histograms), are compared with 
the results relative to the v^v^yy process, proceeding only via Z-boson exchange (dashed 
histograms). For the z/z/77 final state, a cut on the recoil mass around the Z mass (80 GeV 
< M reco u < 120 GeV) has been imposed. For the sake of comparison, the SM predictions 
(dotted histograms) normalized to the same luminosity are also shown. The difference 
between the solid and dashed histograms indicates a very moderate contribution due to 
W-Z interference, even for the large a and a c values considered in the simulation. 

A similar analysis is shown in Fig. 5 for a 4/ + 7 final state, as computed by means of 
the exact calculation of WRAP, in comparison with the WWy approximation considered in 
the literature [11, 13]. The sensitivity to the anomalous couplings a ,a c ,a n is shown for 
the 4/ + 7 final state according to two different event selections: no cuts on the invariant 
masses of the decay products (dotted line) and cuts on the invariant masses of the decay 
products around the W mass, i.e. 75 GeV < M u j^ < 85 GeV, (dashed line), in order 
to disentangle, as much as possible, the contributing Feynman graphs with two final- 
state resonant W bosons. The WWy approximation (solid line) predicts a quite different 
sensitivity with respect to the complete 4/ + 7 calculation. 
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By considering variations of the anomalous couplings within the allowed experimental 
bounds, differences at the ten per cent level are registered between the WW^y approxi- 
mation and the 4/ + 7 prediction, when invariant mass cuts are imposed in the 4/ + 7 
calculation. Notice that, if invariant mass cuts are not considered, the 4/ + 7 cross sec- 
tion grows up by a factor of two with respect to the cross section in the presence of 
cuts. Therefore, the WWj approximation should be employed with the due caution in 
QAGC studies, especially if one takes into account that exact 4/ + 7 generators, such 
as WRAP [18] and RacoonWW [15], incorporate the effects of QAGC and are at disposal for 
such experimental studies. 

For the sake of simplicity, only numerical results in the Born approximation have been 
presented in the paper. Anyway, as already remarked and previously discussed in detail 
in refs. [17, 18], the effects of ISR can be properly simulated by means of NUNUGPV and 
WRAP and should be considered in the experimental analysis since they tend to diminish 
the sensitivity on the QAGC at some per cent level, as proved by explicit numerical 
investigation. 

4 Conclusions 

The search for QAGC in radiative events at LEP demands precise predictions for the 
processes e + e~ — > z/z/77, e + e~ — > ^^77 and e + e~ — > 4 fermions+7. To this end, exact cal- 
culations of such processes have been presented, incorporating the contribution of QAGC 
and the large effect of ISR. 

By means of the exact calculations, comparisons with approximate results existing in 
the literature have been performed. It turns out that, for the z/z/77 an d 9577 final states, 
the Z77 approximation works well, being the interference effects present in the complete 
calculations confined at the per cent level, if appropriate cuts around the Z-boson mass 
are required. It has been also demonstrated that the z/z/77 final state with appropriate 
cuts on the recoil mass can be successfully exploited to extract limits on neutral QAGC, 
as a complementary channel to the qcfi^i one. As far as the 4 fermions+7 final states are 
concerned, significant differences are seen between the exact calculation and the WW^f 
approximation, even in the presence of invariant mass cuts around the VF-boson mass. 

In the spirit of the present study, a phenomenological analysis of anomalous gauge 
couplings in radiative events at the energies of future linear colliders, including the effects 
of beam polarization and beamsstrahlung, will be given elsewhere. 

The exact calculations addressed in the present paper are available in the form of 
improved versions of the event generators NUNUGPV and WRAP 1 , which can be used for full 
experimental simulations. 



1 The Monte Carlo programs NUNUGPV and WRAP can be downloaded from the WEB site 
http : //decuxl . pv . inf n. it/~ nicrosi/programs .html 
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